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Abstract 

The enhancement of neutrino oscillations in matter is briefly reviewed. Exact and 
approximate solutions of the equations describing neutrino oscillations in matter 
are discussed. The role of stochasticity of the media that the neutrinos propagate 
through is elucidated. 
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1 Introduction 


Particle and nuclear physicists devoted an increasingly intensive effort dur¬ 
ing the last few decades to searching for evidence of neutrino mass. Recent 
announcements by the Superkamiokande collaboration of the possible oscilla¬ 
tion of atmospheric neutrinos [1] and very high statistics measurements of the 
solar neutrinos [2] brought us one step closer to understanding the nature of 
neutrino mass and mixings. Experiments imply that neutrino mass is small 
and the seesaw mechanism [3], to the development of which Dick Slansky 

* Expanded version of a talk at the Slansky Memorial Symposium, Los Alamos, 
May 1998. 

^ Electronic address: baha@nucth.physics.wisc.edu 
^ Permanent Address 


Preprint submitted to Elsevier Preprint 


1 February 2008 





contributed, is perhaps the simplest model which leads to a small neutrino 
mass. 

If the neutrinos are massive and different flavors mix they will oscillate as 
they propagate in vacuum [4], Dense matter can signihcantly amplify neutrino 
oscillations due to coherent forward scattering. This behavior is known as the 
Mikheyev-Smirnov-Wolfenstein (MSW) effect [5]. Matter effects may play an 
important role in the solar neutrino problem [6-8]; in transmission of solar 
[8] and atmospheric neutrinos [9,10] through the Earth’s core; and shock re¬ 
heating [11] and r-process nucleosynthesis [12] in core-collapse supernovae. If 
the neutrinos have magnetic moments matter effects may also enhance spin- 
flavor precession of neutrinos [13]. 

The equations of motion for the neutrinos in the MSW problem can be solved 
by direct numerical integration, which must be repeated many times when 
a broad range of mixing parameters are considered. This often is not very 
convenient; consequently various approximations are widely used. Exact or 
approximate analytic results allow a greater understanding of the effects of 
parameter changes. The purpose of this article is to present a review of the so¬ 
lutions of the neutrino propagation equations in matter. Recent experimental 
developments and astrophysical implications of the neutrino mass and mixings 
are beyond the scope of this article. Very rapid developments make a medium 
such as the World Wide Web more suitable for the former and the latter was 
recently reviewed elsewhere [14]. Recent experimental developments can be ac¬ 
cessed through the special home page at SPIRES [15] and theoretical results 
at the Institute for Advanced Study [16] and the University of Pennsylvania 
[17]. An assessment of the Superkamiokande solar neutrino data was recently 
given by Bahcall, Krastev, and Smirnov [6]. A number of recent reviews cover 
implications of recent results for neutrino properties [18]. 


2 Outline of the MSW Effect 


The evolution of flavor eigenstates in matter is governed by the equation [5,20] 
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for the mixing of two active neutrino flavors and 
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for the active-sterile mixing. In these equations 
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Sm'^ = m\ — m\ is the vacuum mass-squared splitting, 9^ is the vacuum mixing 
angle, Gp is the Fermi constant, and Ne{x) and Nn{x) are the number density 
of electrons and neutrons respectively in the medium. 

In a number of cases adiabatic basis greatly simplihes the problem. By making 
the change of basis 
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the flavor-basis Hamiltonian of Eq. (1) can be instantaneously diagonalized. 
The matter mixing angle in Eq. (5) is dehned via 
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In the adiabatic basis the evolution equation takes the form 

IK + (p‘^{x) —ihd'{x) 
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where prime denotes derivative with respect to x. Since the 2x2 “Hamiltonian” 
in Eq. (8) is an element of the SU{2) algebra, the resulting time-evolution 
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operator is an element of the SU{2) group. Hence it can be written in the 
form [19] 




(9) 


where and are solutions of Eq. ( 8 ) with the initial conditions 

= 1 and '^ 2 {xi) = 0 . if the matter mixing angle, 0{x), is changing very 
slowly (i.e., adiabatically) its derivatives in Eq. ( 8 ) can be set to zero. In this 
approximation the “Hamiltonian” in the adiabatic basis is diagonal and the 
system remains in one of the matter eigenstates. 

To calculate the electron neutrino survival probability Eq. (1) needs to be 
solved with the initial conditions Te = 1 and = 0. Using Eq. (9) the 
general solution satisfying these initial conditions can be written as 


'^e{x) =cos 6 *(a;)[cos 6 'jTi(a;) — sin 6 'i\h 2 (a;)] 

+ sin 6 '(a;)[cos 6 *i\h 2 (a;) + sin 6 *jT][(a;)], ( 10 ) 

where 9i is the initial matter angle. Once the neutrinos leave the dense matter 
(e.g. the Sun), the solutions of Eq. ( 8 ) are particularly simple. Inserting these 
into Eq. (10) we obtain the electron neutrino amplitude at a distance L from 
the solar surface to be 


'he(T) = cos 6 '^[cos 6 'i\['i^( 5 ) - sin exp 

+ sin 6 '^[cos 6 'i\I' 2 ,( 5 ) - sin^jT); exp 



( 11 ) 


where Ti (s) and 'h 2 ,(s') are the values of Ti(a;) and 'I' 2 (a;) on the solar surface. 
The electron neutrino survival probability averaged over the detector position, 
L, is then given by 


P(z/e ^ Ve) = {\^e{L)\^)L = ^ + ^ COS 261^ COS 26'* (l - 2|T2,(5)P) 

- - cos 26, sin 26, (Ti,(s)T 2 ,(s) + ^t( 5 )^ 2 ,(S)) • ( 12 ) 

If the initial density is rather large, then cos 20* ~ — 1 and sin 20* ~ 0 and the 
last term in Eq. (12) is very small. Different neutrinos arriving the detector 
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carry different phases if they are produced over an extended source. Even if 
the initial matter density is not very large, averaging over the source position 
makes the last term very small as these phases average to zero. The completely 
averaged result for the electron neutrino survival probability is then given by 
[ 21 ] 

1 1 

P{j^e ^ ^e) 2 2 2/\op) 1 (12) 

where the hopping probability is 

Tiiop = |tE'2,(S)|^, (14) 

obtained by solving Eq. (8) with the initial conditions Ti(a;i) = 1 and ^ 2 {,Xi) = 
0. Note that, since in the adiabatic limit T 2 ,(s’) remains to be zero Phop = 0. 


3 Exact Solutions 


Exact solutions for the neutrino propagation equations in matter exist for a 
limited class of density prohles that satisfy an integrability condition called 
shape invariance [22]. To illustrate this integrability condition we introduce 
the operators 
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Using Eq. (15) Eq. (1) takes the form 


i_Te(a;) = yAT^(a;), 

i+T^(a;) = VATe(a;). ( 16 ) 

The shape invariance condition can be expressed in terms of the operators 
defined in Eq. (15) [23] 

yi_(ai)/i-i_(ai) = /i-|_(a2)^-(®2) + P(ai). (11") 


We also introduce a similarity transformation which formally replaces Oi by 
0.2 ■ 

f{a,)0{a,)f-\a,) = ^( 02 ). (18) 
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The MSW equations take a particularly simple form using the operators [24] 


B+ = A+{ai)f{ai) 

B_ =T~^{ai)A_{ai) , (19) 

which satisfy the commutation relation: 

[B_,B+]=R{ao), ( 20 ) 

where oq is dehned using the identity 

/2(an) = T(ai)i?(a„_i)T"^(ai), (21) 

with n = 1. Two additional commutation relations 

[B^B_, Bl] = (i?(ai) + R{a2) + ■ ■ +i?(a„))5”, (22) 

and 

[B+B_,BzA = (^(ai) + R{a2) + ■ ■ +i?(a„))5Z" , (23) 


can easily be proven by induction. 

Using the operators introduced in Eq. (19), Eq. (1) can be rewritten as 

B+B_^,{x) = ATe(a;). (24) 

Eqs. (22) and (23) suggest that 5+ and B_ can be used as ladder operators 
to solve Eq. (24). Introducing 

T® ~ exp i j Lp{x; ai)dx^. (25) 

one observes that 

i_(ai)T® = 0 = B_¥°\ (26) 

If the function 
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can be analytically continued so that the condition 


/(fi) = A 


( 28 ) 


is satished for a particular (in general, complex) value of fi, then Eq. (22) 
implies that one solution of Eq. (24) is Similarly the wavefunction 




( 0 ) 


exp J (p{x; ao)dx 


(29) 


satishes the equation 

= 0 . 


(30) 


Then a second solution of Eq. (24) is given by 51^”^Hence for shape in¬ 
variant electron densities the exact electron neutrino amplitude can be written 
as [24] 


'he(a^) = /3.B+ exp i J ip{x]ai)dx^ 

+'jBZ^~^ exp (^+i J Lp{x; ao)dx^, (31) 

where (3 and 7 are to be determined using the initial conditions 'hi(a;j) = 1 
and = 0 . 

For the linear density prohle 

Ne{x) = No - Nq{x - xr), (32) 

where Nq is the resonant density: 

2 V 2 GfNqE = 5m^ cos29^, (33) 


using the technique described above we can easily write down the hopping 
probability 

-Phop = \^ 2 {xf)\‘^ = exp (-ttH) , (34) 


where 


5m^ sin^ 29^ Nq 
AE cos 20^ iVg 


( 35 ) 
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This is the standard Landau-Zener result [21,25]. 
For the exponential density prohle 




( 36 ) 


where Nq is the resonant density given in Eq. (33), the hopping probability is 

|26] 
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where we dehned 
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4 Supersymmetric Uniform Approximation 


The coupled hrst-order equations for the flavor-basis wave functions can be 
decoupled to yield a second order equation for only the electron neutrino 
propagation 


h- 






dx^ 


A + (p^(x) + ih(p'(x) ^e(x) = 0. 


(39) 


The large body of literature on the second-order differential equations of math¬ 
ematical physics motivates using a semiclassical approximation for the solu¬ 
tions of Eq. (39). The standard semiclassical approximation gives the adiabatic 
evolution [27]. For a monotonically changing density prohle supersymmetric 
uniform approximation yields [28] 


Fhop = exp(-7rQ), 
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i 5m^ 
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where Tq and tq are the turning points (zeros) of the integrand. In this expres¬ 
sion we introduced the scaled density 
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Fig. 1. The electron neutrino survival probability for the Sun [28]. The solid line 
is calculated using Eq. (40). The dashed line is the exact (numerical) result. The 
dotted line is the linear Landau-Zener result. In the top figure, the lines are indis¬ 
tinguishable. An exponential density with parameters chosen to approximate the 
Sun was used [29]. 


where Ne is the number density of electrons in the medium. By analytic con¬ 
tinuation, this complex integral is primarily sensitive to densities near the 
resonance point. The validity of this approximate expression is illustrated in 
Figure 1. As this hgure illustrates the approximation breaks down in the ex¬ 
treme non-adiabatic limit (i.e., as 5w? 0). Hence it is referred to as the 

quasi-adiabatic approximation. 


The near-exponential form of the density prohle in the Sun [29] motivates an 
expansion of the electron number density scale height, r^, in powers of density: 


- Ts 


N,(r) 

Kir) 


T.f’nK, 


(42) 


where prime denotes derivative with respect to r. In this expression a mi¬ 
nus sign is introduced because we assumed that density prohle decreases a 
r increases. (For an exponential density prohle, N^. ~ only the n = 0 

term is present). To help assess the appropriateness of such an expansion the 
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Fig. 2. Electron number density scale height (cf. Eq. (42)) as a function of the 
radius for the Sun [29]. The dashed line is the exponential fit over the whole Sun. 
The shaded are indicates where the small angle MSW resonance takes place for 
neutrinos with energies 5 < E < 15 MeV. 


density scale height for the Sun calculated using the Standard Solar Model 
density prohle is plotted in Figure 2. One observes that there is a signihcant 
deviation from a simple exponential prohle over the entire Sun. However the 
expansion of Eq. (42) needs to hold only in the MSW resonance region, indi¬ 
cated by the shaded area in the hgure. Real-time counting detectors such as 
Superkamiokande and Sudbury Neutrino Observatory, which can get informa¬ 
tion about energy spectra, are sensitive to neutrinos with energies greater than 
about 5 MeV. For the small angle solution (sin 26* ~ 0.01 and = 5 x 10“® 
eV^), the resonance for a 5 MeV neutrino occurs at about 0.35 Rq and for a 
15 MeV neutrino at about 0.45 Rq (the shaded area in the hgure). In that 
region the density prohle is approximately exponential and one expects that 
it should be sufficient to keep only a few terms in the expansion in Eq. (42) 
to represent the density prohle of the Standard Solar Model. 

Inserting the expansion of Eq. (42) into Eq. (40), and using an integral repre¬ 
sentation of the Legendre functions, one obtains [30] 


SlTl^ ( 

Q = —— |&o(l - cos26*^) 

{2V2Gfe) 2^ [^n-i(cos20„) - P„+i(cos20,)]} , (43) 
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(MeV) 

Fig. 3. (a) Spectrum distortion at SNO for the small-angle MSW solution 
{5rn? ~ 5 X 10“® eV^ and sin20 ~ 0.01). The solid line is the exact numerical 
solution. The dashed, dot-dashed, and dotted lines result from values of n up to 0, 
1, and 2 in Eq. (43). The error bars on the exact numerical result correspond to 
two and five years of data collection. The dot-dot-dot-dashed line is the spectrum 
without MSW oscillations, normalized to the same total rate as with MSW oscil¬ 
lations. Note that on the scale of this figure the n = 1 and n = 2 lines are not 
distinguishable from the exact answer, (b) The relative error arising from the use 
of Eq. (43). 


where Pn is the Legendre polynomial of order n. The n = 0 term in Eq. (43) 
represents the contribution of the exponential density profile alone. Eq. (43) 
directly connects an expansion of the logarithm of the hopping probability in 
powers of 1/E to an expansion of the density scale height. That is, it provides 
a direct connection between iVe(r) and Pu{Ey). Eq. (43) provides a quick and 
accurate alternative to numerical integration of the MSW equation for any 
monotonically-changing density profile for a wide range of mixing parameters. 
The accuracy of the expansion of Eq. (43) is illustrated in Figure 3 where 
the spectrum distortion for the small angle MSW solution is plotted. In this 
calculation we used the method of Ref. [31] and neglected backgrounds. The 
neutrino-deuterium charged-current cross-sections were calculated using the 
code of Bahcall and Lisi [32]. One observes that for the Sun, where the density 
prohle is nearly exponential in the MSW resonance region, the hrst two terms 
in the expansion provide an excellent approximation to the neutrino survival 
probability. 
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5 Neutrino Propagation in Stochastic Media 


In implementing the MSW solution to the solar neutrino problem one typi¬ 
cally assumes that the electron density of the Sun is a monotonically decreas¬ 
ing function of the distance from the core and ignores potentially de-cohering 
effects [33]. To understand such effects one possibility is to study paramet¬ 
ric changes in the density or the role of matter currents [34], In this regard, 
Loreti and Balantekin [35] considered neutrino propagation in stochastic me¬ 
dia. They studied the situation where the electron density in the medium has 
two components, one average component given by the Standard Solar Model or 
Supernova Model, etc. and one fluctuating component. Then the Hamiltonian 
in Eq. (1) takes the form 

H = cos 20 + ^G^Neir) + iV:(r))j a. -F sin20j a..(44) 


where one imposes for consistency 


mr)} = 0 , 


(45) 


and a two-body correlation function 


{N'(r)N'(r')) = 0‘ NJr) ]V,(r') exp(-|r - r'l/Xc). 


(46) 


In the calculations of the Wisconsin group the fluctuations are typically taken 
to be subject to colored noise, i.e. higher order correlations 

/i 2 ... = (iv:(n)iv:(r 2 )---) (47) 


are taken to be 


/i234 — / 12/34 + / 13/24 + / 14 / 23 ) 


(48) 


and so on. 

Mean survival probability for the electron neutrino in the Sun is shown in 
Figure 4 [36] where fluctuations are imposed on the average solar electron 
density given by the Bahcall-Pinsonneault model. 

One notes that for very large fluctuations complete flavor de-polarization 
should be achieved, i.e. the neutrino survival probability is 0.5, the same as the 
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Fig. 4. Mean electron neutrino survival probability in the sun with fluctuations. 
The average electron density is given by the Standard Solar Model of Bahcall and 
Pinsonneault [37] and sin^ 26 = 0.01. 


vacuum oscillation probability for long distances. To illustrate this behavior 
the results from the physically unrealistic case of 50% fluctuations are shown. 
Also the effect of the fluctuations is largest when the neutrino propagation in 
their absence is adiabatic. This scenario was applied to the neutrino convection 
in a core-collapse supernova where the adiabaticity condition is satisfied [38]. 
Similar results were also obtained by other authors [39-42]. It may be possible 
to test solar matter density fluctuations at the BOREXINO detector currently 
under construction [43]. Propagation of a neutrino with a magnetic moment 
in a random magnetic moment has also been investigated [35,44]. Also if the 
magnetic held in a polarized medium has a domain structure with different 
strength and direction in different domains, the modihcation of the potential 
felt by the neutrinos due polarized electrons will have a random character [46]. 
Using the formalism sketched above, it is possible to calculate not only the 
mean survival probability, but also the variance, a, of the huctuations to get 
a feeling for the distribution of the survival probabilities [36] as illustrated in 
Figure 5. 


In these calculations the correlation length r is taken to be very small, of the 
order of 10 km., to be consistent with the helioseismic observations of the 
sound speed [45]. In the opposite limit of very large correlation lengths are 
very interesting result is obtained [38], namely the averaged density matrix is 


13 






Fig. 5. Mean electron neutrino survival probability plus minus a in the sun with 
fluctuations. The average electron density is given by the Standard Solar Model of 
Bahcall and Pinsonneault and sin^ 20 = 0.01. Panels (a), (b), (c), and (d) correspond 
to an average fluctuation of 1%, 2%, 4%, and 8% respectively. 


given as an integral 


lim (p(r)) 

Tc^OO 



dx exp [—x‘^/ (2/9^)]/5(r, x), 


(49) 


reminiscent of the channel-coupling problem in nuclear physics [47]. Even 
though this limit is not appropriate to the solar fluctuations it may be ap¬ 
plicable to a number of other astrophysical situations. 
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